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This  phenomenon  of  internal  rupture  is  closely  related  to  the  veil -known 
phenomenon  of  cavitation  in  fluids,  and  has  therefore  been  dubbed  "cavitation" 
even  in  solids.  Since  the  nucleation  of  a  void  in  a  stressed  solid  is  often  a 
precursor  to  failure,  it  would  be  useful  to  be  able  to  predict  the  conditions 
under  which  cavitation  occurs. 

The  present  project  has  been  aimed  at  trying  to  determine  whether  cavita¬ 
tion,  viewed  as  an  intrinsic  material  instability,  provides  insight  into  the 
phenomenon  of  void  nucleation.  The  project  had  two  principal  goals:  one,  to  bet¬ 
ter  understand  the  phenomenon  of  spherically  (of  cylindrically)  symmetric  cavi¬ 
tation,  and  two,  to  study  the  phenomenon  of  cavitation  in  circumstances  that  do 
not  possess  such  symmetry.  This  report  contains  the  results  of  the  investigation. 
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DESCRIPTION  OF  PROBLEM  AND  RESULTS : 

Consider  a  solid  sphere  of  unit  radius  subjected  to  a  uniformly- distributed 
radial  tensile  traction  p  on  its  surface.  In  response  to  this  loading,  the  body 
can  expand  into  a  larger  solid  sphere.  However  for  certain  materials,  when  p  is 
sufficiently  large,  it  is  energetically  more  favorable  for  the  body  to  instead 
develop  an  internal  spherical  cavity.  A  bifurcation  analysis  yields  the  value 
pcr  of  the  applied  stress  at  which  this  instability  occurs,  without  the  need  for 
an  ad  hoc  failure  criterion,  Ball(1982). 

This  phenomenon  of  internal  rupture  is  closely  related  to  the  well-known 
phenomenon  of  cavitation  in  fluids,  and  has  therefore  been  dubbed  "cavitation" 
even  in  solids.  Since  the  nucleation  of  a  void  in  a  stressed  solid  is  often  a 
precursor  to  failure,  it  would  be  useful  to  be  able  to  predict  the  conditions 
under  which  cavitation  occurs . 

The  present  project  has  been  aimed  at  trying  to  determine  whether  cavita¬ 
tion,  viewed  as  an  intrinsic  material  instability,  provides  insight  into  the 
phenomenon  of  void  nucleation.  The  project  had  two  principal  goals:  one,  to  bet¬ 
ter  understand  the  phenomenon  of  spherically  (of  cylindrically)  symmetric  cavi¬ 
tation,  and  two,  to  study  the  phenomenon  of  cavitation  in  circumstances  that  do 
not  possess  such  symmetry. 

PART  A:  RADIALLY  SYMMETRIC  CAVITATION 

The  effect  of  compressibility  and. hardening  on  cavitation: 

The  study  of  cavitation  in  compressible  materials  is  considerably  more  dif¬ 
ficult  than  the  corresponding  problem  for  incompressible  materials  because  of 
the  difficulty  involved  in  solving  the  associated  nonlinear  differential  equa¬ 
tion.  We  studied  this  problem  for  a  special  class  of  compressible  elastic  mate¬ 
rials,  the  material  model  being  a  generalization  of  one  for  foam  rubber,  and  the 
aim  being  to  examine  the  qualitative  influence  of  material  compressibility  on 
cavitation.  Specifically,  a  bifurcation  problem  for  a  solid  circular  cylinder 
composed  of  this  elastic  material  was  studied.  The  curved  surface  of  the  cylin¬ 
der  was  subjected  to  a  radial  stretch  A  (>1) ,  and  the  cylinder  was  in  a  state  of 
plane  strain.  A  cylindrical  cavity,  coaxial  with  the  cylinder,  is  found  to 
emerge  when  A  reaches  a  critical  value  (say  *cr)  at  which  the  homogeneous 
deformation  becomes  unstable.  We  find  that  Acr  decreases  as  the  material  becomes 
softer  and  also  as  it  becomes  less  compressible.  The  corresponding  value  of 
radial  true  stress  rcr  also  decreases  as  the  material  becomes  softer,  but 
increases  as  it  becomes  less  compressible.  Figures  1  and  2  show  this  variation. 

Reference:  N.  Ertan,  Influence  of  compressibility  and  hardening  on  cavitation, 
ASCE  Journal  of  Engineering  Mechanics,  Vol.  114,  (1988),  pp.  1231-1244. 

The  effect  of  rate -dependence  on  the  growth  of  an  infinitesimal  cavitv: 

In  this  study  we  examined  the  effect  of  rate -dependence  on  the  phenomenon  of 
cavitation  in  an  incompressible  material.  Specifically,  a  sphere  containing  a 
traction-free  void  of  infinitesimal  initial  radius  was  considered,  and  the  outer 
surface  of  the  sphere  was  subjected  to  a  prescribed  uniform  radial  nominal 
stress  p,  which  was  suddenly  applied  and  then  held  constant.  The  sphere  was  com¬ 
posed  of  a  particular  class  of  incompressible  rate -dependent  materials.  The 
large  strains  which  occur  in  the  vicinity  of  the  void  were  accounted  for  in  the 
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analysis  and  the  problem  was  reduced  to  a  nonlinear  initial -value  problem,  which 
was  then  studied  qualitatively  through  a  phase-plane  analysis.  The  principal 
result  derived  in  this  study  is  characterized  by  two  equations  relating  the 
applied  stress  p  and  the  current  cavity  radius  b:  p  -  p(b)  and  p  -  p(b) .  The 
first  of  these  describes  a  curve  that  separates  the  (p,b) -plane  into  regions 
where  cavitation  does  and  does  not  occur.  The  second  describes  a  curve  which 
subdivides  the  former  subregion  --  the  post-cavitation  region  --  into  domains 
where  void  growth  occurs  slowly  and  rapidly.  Figure  3  shows  an  example  of  these 
curves  for  a  particular  set  of  material  parameters. 

Reference:  R.  Abeyaratne  and  H-s.  Hou,  Growth  of  an  Infinitesimal  Cavity  in  a 
Rate -Dependent  Solid,  ASME  Journal  of  Applied  Mechanics,  Vol.  56,  (1989),  pp. 
40-46. 


The  effect  of  axial  stretch  on  cavitation  in  an  elastic  cylinder 

In  this  study  the  phenomenon  of  cavitation  in  an  elastic  cylinder  subjected 
to  a  combined  axial  stretch  Az  and  a  radial  traction  p  was  examined.  One 
possible  response  of  the  cylinder  to  this  loading,  for  all  Az  and  p,  is  a  pure 
homogeneous  deformation.  However,  for  some  materials,  the  homogeneous  deforma¬ 
tion  becomes  unstable  at  certain  critical  values  of  the  pair  (Az,p)  and  a  second 
deformation  involving  a  cylindrical  cavity  emerges.  We  determined  the  values  of 
(Az,p)  at  which  this  happens.  The  results  are  displayed  by  constructing  a  curve 
in  the  (Az,p) -plane  which  divides  this  plane  into  regions  where  the  homogeneous 
deformation  is  stable  and  unstable. 

From  the  results  of  this  study  we  have  observed  that  a  material  which  does 
not  exhibit  cavitation  at  one  value  of  the  axial  stretch  Az  might  do  so  when  it 
is  subjected  to  some  other  value  of  Az.  We  also  observed  that  a  material  which, 
at  some  Az,  possesses  cavitated  deformations  which  are  unstable,  could  exhibit 
stable  cavitated  deformations  when  a  different  axial  stretch  is  imposed.  Figure 
4  shows  the  critical  curves  in  the  (Az,p) -plane  for  different  values  of  material 
hardening. 

Reference:  H-s.  Hou  and  Y.  Zhang,  The  effect  of  axial  stretch  on  cavitation, 
accepted  for  publication  in  the  International  Journal  of  Nonlinear  Mechanics, 
November  1989. 


Comparison  of  analysis  with  experiments 

Gent  and  Lindley(1958)  conducted  experiments  on  thin  disks  of  vulcanized 
rubber  by  bonding  its  faces  to  two  stiff  metal  disks,  and  then  applying  uni¬ 
axial  loads  to  the  metal  disks  normal  to  its  faces.  They  observed  direct  and 
indirect  evidence  of  flaw  initiation  at  the  center  of  the  rubber  specimen  when 
the  load  reached  a  critical  value.  By  repeating  the  experiment  on  a  series  of 
specimens  which  had  been  pre-treated  differently,  they  were  able  to  relate  the 
critical  load  at  void  initiation  to  the  Young's  modulus  of  the  material. 

Oberth  and  Bruenner(1965)  conducted  tensile  tests  on  a  standard  specimen  of 
polyurethane  in  which  they  had  embedded  a  small  metallic  sphere;  the  sphere  was 
well  bonded  to  the  surrounding  matrix  material.  During  the  tension  tests  they 
observed  the  appearance  of  small  cavities  (near  the  poles  of  the  spherical  par- 
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ticle)  when  the  load  reached  a  critical  value.  They  too  repeated  the  tests  on  a 
series  of  specimens  that  had  been  pre- treated  differently,  and  were  able  to 
relate  the  Young's  modulus  of  the  material  to  the  critical  load  at  void  initia¬ 
tion  . 

We  carried  out  a  finite  element  analysis  of  the  two  test  specimens  used  in 
these  experiments,  and  used  the  results  of  these  numerical  studies,  in  conjunc¬ 
tion  with  the  predictions  of  a  cavitation  analysis,  to  predict  theoretically  the 
load  when  cavitation  occurs;  in  both  cases,  the  results  were  in  striking  agree¬ 
ment  with  the  experimental  observations. 

Reference:  R.  Stringfellow  and  R.  Abeyaratne,  Cavitation  in  an  elastomer:  com¬ 
parison  of  theory  with  experiment,  Journal  of  Materials  Science  and  Engineering, 
Vol .  A112,  (1989),  pp.  127-131. 


Void  collapse 

In  this  study  we  examined  the  collapse  of  a  void,  the  phenomenon  in  which  an 
existing  void  closes  due  to  an  instability.  This  is  the  opposite  phenomenon  to 
cavitation.  A  complete  analysis  of  the  collapse  of  a  void  in  an  incompressible 
elastic  material  was  carried  out  in  the  special  cases  of  spherical  and  axial 
symmetry  .  The  analysis  predicts  the  critical  load  at  which  collapse  occurs. 

We  considered  large,  radially  symmetric,  deformations  of  an  infinite  medium 
composed  of  an  isotropic,  incompressible  elastic  material  surrounding  a  trac¬ 
tion-free  spherical  cavity.  The  body  is  subjected  to  a  uniform  pressure  at  infi¬ 
nity  and  we  examine  the  possibility  of  void  collapse,  i.e.  the  possibility  that 
the  void  radius  becomes  zero  at  a  finite  value  of  the  applied  stress.  This  does 
not  occur  in  all  materials.  We  have  determined  the  complete  sub-class  of  incom¬ 
pressible  elastic  materials  that  do  exhibit  this  phenomenon,  and  for  such  mate¬ 
rials,  found  the  value  of  the  applied  load  at  which  collapse  occurs.  We  have 
also  examined  the  stability  of  the  deformation,  both  from  minimum  potential 
energy  and  dynamic  stability  points  of  view.  The  results  were  illustrated  by 
considering  a  particular  class  of  power- law  materials.  Figure  3  shows  the  varia¬ 
tion  of  the  cavity  radius  b  with  the  applied  pressure  p  for  three  values  of  the 
hardening  exponent. 

In  certain  respects,  the  results  for  void  collapse  are  complementary  to 
those  for  void  expansion.  In  this  complementary  problem,  one  is  interested  in 
examining  if,  and  when,  a  cavity  In  an  infinite  medium  can  grow  without  bound. 
Many  formal  similarities  exist  between  the  results  of  the  analyses  of  these  two 
problems . 

Reference:  R.  Abeyaratne  and  H-s.  Hou,  Void  collapse  in  an  incompressible  elas¬ 
tic  solid,  accepted  for  publication  in  the  Journal  of  Elasticity.  July  1989. 


The  relationship  between  local  and  global  instabilities  during  cavitation. 

Consider  a  solid  sphere  subjected  to  a  uniform  radial  tensile  load,  which, 
at  some  critical  value  of  the  load,  bifurcates  into  a  configuration  involving  a 
concentric  internal  cavity.  Consider  an  infinitesimal  cube  of  material  within 
this  body,  just  prior  to  cavitation,  which  is  oriented  with  its  edges  parallel 
to  the  radial  and  hoop  directions.  In  this  configuration,  this  material  element 
is  subjected  to  a  deformation  with  equal  principal  stretches  “  *2  ”  *3  anci  a 
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hydrostatic  state  of  stress.  Immediately  after  cavitation,  this  same  material 
element  has  >.ne  geometric  shape  of  a  plate  with  A^  <  1  and  *2  ”  *3  (where  A^ 
is  the  radial  stretch  and  A2,  A3  are  the  hoop  stretches);  see  Figure  6.  Thus  it 
is  apparent  that  when  the  body  undergoes  the  (global)  cavitation  instability, 
each  infinitesimal  material  el  nt  undergoes  this  local  "cube  -*  p_ate"  insta¬ 
bility  (which  we  refer  to  as  an  asymmetric  instability) .  This  suggests  that 
cavitation  might  be  viewed  as  the  occurrence  of  a  continuum  of  local  asymmetric 
instabilities,  so  that 

i)  a  material  that  exhibits  the  cavitation  instability  must  necessarily  exhibit 
the  asymmetric  instability; 

ii)  the  critical  loan  at  cavitation  must  be  a  suitable  average  of  a  sequence  of 
critical  loads  for  the  asymmetric  instability^ 

iii) a  material  that  exhibits  stable  cavitation  must  necessarily  exhibit  stable 
asymmetric  deformations. 

We  have  been  able  to  prove  these  results,  thereby  allowing  us  to  understand  the 
intrinsic  nature  of  the  cavitation  instability  in  terms  of  a  more  simple  insta¬ 
bility  that  occurs  locally. 

Reference:  H-s.  Hou  and  R.  Abeyaratne ,  On  the  cavitated  instability  and  its 
relationship  to  a  continuum  of  local  asymmetric  conf iguracions? ,  manuscript  in 
preparation . 


PART  B:  NON -RADIALLY  SYMMETRIC  CAVITATION 

Most  previous  studies  on  cavitation,  including  those  described  in  Part  A 
above,  have  been  restricted  to  the  case  of  radial  (i.e.  spherical  or  cylindri¬ 
cal)  symmetry.  Many  practical  situations  do  not  possess  such  symmetry  and  it  is 
therefore  important  to  extend  the  radially- symmetric  studies  to  more  general 
situations.  The  particular  three-dimensional  case  of  axial - symmetr; ,  where  the 
nucleated  void  is,  roughly,  ellipsoidal  in  shape,  is  a  special  case  of  some 
importance . 

Consider  an  incompressible,  isotropic  medium,  subjected  to  a  remotely 
applied  axially- symmetric  state  of  stress;  let  denote  the  axial  stress  and  r 2 
each  of  the  two  equal  in-plane  stresses.  The  body  deforms  in  an  axially - 
symmetric  manner.  If  -  r^,  we  have  spherical  symmetry  and  an  (infinitesimal) 
spherical  void  appears  at  the  origin,  when  the  applied  stress  reaches  the  value 
rl  “  r2  ”  rcr  t^ie  general  case  (r^  h  t 2)  one  expects  to  have  a  function 
C r 1 ,  r2)  characterizing  a  failure  or  "cavitation  curve"  such  that  when  ^(r^, 

<  0  the  body  deforms  homogeneously,  when  r 2)  -  0  an  infinitesimal 

(not-necessarily-spherical)  void  is  nucleated,  and  when  ^(r^,  r 2)  >  0  the  cavity 
expands;  in  order  to  conform  with  the  symmetric  case,  one  must  have  <^(rcr,  rcr ) 

-  0.  Determining  this  cavitation  curve  is  the  focus  of  the  studies  described  in 
this  part. 

First,  we  showed  that  the  phenomenon  of  cavitation  is  always  preceded  by  a 
different  type  of  instability  if  the  prescribed  loading  is  dead,  and  therefore, 
in  all  subsequent  work  we  restricted  attention  on  the  case  of  prescribed  dis¬ 
placement  or  prescribed  true  stress.  Next,  we  determined  an  exact  analytical 
expression  for  the  cavitation  function  ^  for  a  special  compressible  elastic 
material.  We  then  developed  an  analytical  procedure  which  led  to  an  approximate 
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Reference:  R.  Abeyaratne  and  H-s.  Hou,  On  the  occurrence  of  the  cavitation 
instability  relative  to  the  asymmetric  instability  under  symmetric  dead  load 
conditions,  submitted  to  Quarterly  Journal  of  Mechanics  and  Applied  Mathematics, 
April  1990. 


Cavitation  and  void  collapse  under  plane  strain  bi-axial  loading  in  a  special 
compressible  elastic  material: 

Consider  a  state  of  plane  strain  of  an  infinite  medium  surrounding  a  trac¬ 
tion-free  circular  hole,  and  suppose  that  the  medium  is  composed  of  the  power- 
law  material  considered  by  Varley  and  Cumberbatch(1980) :  in  such  a  material,  the 
nominal  stress  a  and  the  axial  stretch  A  are  related,  in  plane  strain  uni-axial 
tension,  by 

o  -  (E/k)  (A*'1  -  A'1}, 

where  E  and  k  are  material  constants.  At  infinity,  the  body  is  subjected  to  a 
state  of  bi-axial  stress 

-*  ,  a22  a2>  a12  an<^  a2\ 

We  have  analyzed  the  associated  problem  and  determined  the  resulting  deformation 
and  stress  fields.  Suppose  for  example  that  the  hardening  exponent  has  the 
value  k-2.  Figure  7  shows  a  sketch  of  the  ,02) -plane .  The  two  bold  curves  in 

the  figure  form  a  critical  boundary.  If  the  point  corresponding  to  the 

remotely  applied  stress  lies  in  the  hatched  region,  the  cavity  is  open  in  the 
deformed  configuration  and  has  an  elliptical  shape.  As  the  point  (01,02) 
approaches  one  of  the  boundaries  of  the  hatched  region,  the  cavity  collapses 
into  a  crack  as  shown. 

Suppose  instead  that  the  hardening  exponent  in  the  constitutive  law  has  the 
value  k-1.  In  this  case  the  cavity  is  generally  ' 'peanut -shaped '' in  the  deformed 
configuration  as  shown  in  Figure  8.  As  (0^,02)  approaches  either  of  the  two  bold 
lines  with  positive  slope,  the  cavity  collapses  partially  by  contacting  along 
part  of  its  boundary  as  shown  in  the  figure.  On  the  other  hand  when  (0^,02) 
approaches  the  bold  line  with  negative  slope,  the  deformed  cavity  size  becomes 
infinite  corresponding  to  cavitation. 

Reference:  R.  Abeyaratne  and  H-s.  Hou,  Void  collapse  in  an  elastic  solid, 
accepted  for  publication  in  the  Journal  of  Elasticity.  July  1989. 


Analytical  approximation  to  the  cavitation  curve: 

Consider  an  infinite  incompressible  medium,  subjected  at  infinity  to  an  axi- 
symmetric  stress  state  (r^,  r2,  r 2 ).  We  wish  to  determine  the  resulting 
deformation  of  this  body.  One  deformation  that  the  body  can  undergo,  no  matter 
what  the  values  of  and  r2*  is  a  pure  homogeneous  one.  Determining  other 
possible  deformations,  in  exact  analytical  closed  form,  is  a  difficult  task; 
consequently,  in  this  part  of  the  present  study,  we  attempted  to  find  an  approx¬ 
imate  analytical  deformation  field. 

Consider  the  sub-class,  say  ’A.',  of  all  kinematically  admissible  deforma- 
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tions  which  are  of  the  form 
Yl  -  fCrJxi, 

y2  “  S(r)x2,  •  (1) 

y3  -  g(r)x3; 

here  (x3,x2,x3)  are  the  coordinates  of  a  particle  before  deformation  and 
(yi.y2-y3>  are  its  coordinates  in  the  deformed  configuration,  r2  -  xi2+x22+x32, 
and  f  and  g  are  positive-valued  functions.  Note  that  deformations  in  this  class 
possess  axi- symmetry;  moreover,  the  special  case  f-g  corresponds  to  radially 
symmetric  deformations,  while  the  special  case  f-constant,  g-constant  describes 
a  homogeneous  axi- symmetric  deformation.  In  view  of  the  incompressibility  of  the 
material  one  can  show  that  deformations  in  *A'  must  necessarily  have  the  form 

yi  -  A'2  (1  +  03/R3)1/3  x1( 

y2  -  A  {1  +  ^3/R3)1/3  x2,  ►  (2) 

y3  -  A  [1  +  03/R3)V3  x3,  . 

where  A  >  0  and  0  >  0  are  kinematically  undetermined  constants.  When  0-0  this 

describes  a  pure  homogeneous  deformation,  while  when  A-l  it  corresponds  to  a 
radial  deformation  involving  a  cavity  of  current  radius  0.  In  general,  such  a 
deformation  carries  the  origin  in  the  undeformed  configuration  into  an  ellip¬ 
soidal  cavity;  0  is  a  measure  of  the  size  of  the  cavity  while  A  describes  its 
departure  from  spherical  symmetry. 

Thus,  the  collection  of  deformations  ’A'  corresponds  to  a  two-parameter  fam¬ 
ily  of  deformations  and  the  "best"  values  for  the  parameters  A,  0  may  be  found 
by  minimizing  the  appropriate  total  energy  of  the  body.  Such  a  calculation  leads 
to  the  following  two  equations  relating  the  stresses  and  r2  to  the  geometric 
parameters  A  and  0: 

r3  +  2 r2  -  3/(4jt02)  3E/30, 

T2  '  T1  "  3A/(8w(1+03))  3E/3A ; 

here  E(0,A)  is  the  total  strain  energy  stored  in  the  body  corresponding  to  the 
deformation  (2).  At  the  instant  of  cavitation  one  has  0-0;  setting  0-0  in  (3) 
and  then  eliminating  (in  principle)  the  parameter  A  between  those  two  equations, 
leads  to  a  relationship  between  rj.  and  r2  which  is  the  approximate  equation  of 
the  cavitation  curve  ^(r^,r2)-0.  The  curve  in  Figure  9  shows  the  result  for  the 
neo-Hookean  material.  This  same  procedure  can  be  carried  out  for  elastic-plastic 
materials  as  well,  provided  the  principle  of  virtual  work  is  used  in  place  of 
extremizing  the  total  energy.  The  curves  in  Figure  10  show  the  result  for  a 
piecewise  power- law  material  for  different  values  of  the  hardening  exponent. 

Reference:  H-s.  Hou  and  R.  Abeyaratne,  Axi-symmetric  cavitation  in  solids,  manu¬ 
script  in  preparation. 


-  - 14-  - 


Numerical  verification  of  accuracy  of  approximate  cavitation  curves: 

In  order  to  investigate  the  accuracy  of  the  preceding  analytical  approxima¬ 
tion,  we  undertook  a  numerical  study.  By  using  the  large  strain  analysis 
capabilities  of  the  finite  element  program  ABAQUS ,  we  have  carried  out  numerical 
simulations  for  both  the  neo-Hookean  material  and  piecewise  power-law  elastic- 
plastic  materials.  For  computational  purposes,  we  must,  of  course,  deal  with 
finite  sized  regions  and  finite  sized  cavities.  Thus,  we  have  taken  a  "large" 
body  containing  a  "small"  cavity  and  incremented  the  stress  in  steps  (along 
radial  lines  in  stress  space)  and  solved  for  the  various  fields  including  the 
cavity  dimensions.  In  a  graph  of  cavity  "size"  (-  one  half  the  sum  of  major  and 
minor  axes)  versus  stress,  the  cavity  size  is  found  to  growth  very  slowly  first 
and  then  to  suddenly  grow  rapidly.  We  determined  the  limit  load,  i.e.  the  maxi¬ 
mum  value  of  applied  stress,  and  took  it  to  be  an  approxii  ion  for  the  stress 
level  at  cavitation. 

The  circles  in  Figure  9  show  the  cavitation  curve  for  the  neo-Hookean  mate¬ 
rial  as  determined  by  the  finite  element  solution.  We  see  that  the  analytical 
solution  is  a  good  approximation  and  it  certainly  shows  the  correct  trends. 
Figure  11  shows  the  analytical  and  numerical  cavitation  curves  (corresponding  to 
the  curve  and  the  circles  respectively)  for  an  elastic -perfectly  plastic  mate¬ 
rial;  the  circles  in  that  figure  were  taken  from  the  numerical  solution  of  Hut¬ 
chinson  et  al.(1989).  Similar  comparisons  for  materials  with  strain  hardening 
were  also  carried  out  but  are  not  shown  here.  Finally  Figure  12  shows  a  compari¬ 
son  between  the  actual  deformation  fields;  the  two  curves  shown  are  the  deformed 
cavity  shape  according  to  (2)  at  two  different  values  of  load,  while  the  circles 
correspond  to  the  finite  element  solution;  the  figure  is  drawn  for  the  neo- 
Hookean  material. 

Reference:  H-s.  Hou  and  R.  Abeyaratne,  Axi-symmetric  cavitation  in  solids,  manu¬ 
script  in  preparation. 


Application  of  cavitation  curves:  comparison  of  theory  with  experiment. 

We  have  used  the  failure  curves  ^  generated  by  the  aforementioned  studies, 
and  applied  it  to  two  experimental  situations:  one,  the  experiments  of  Gent  and 
Lindley(1958)  on  thin  rubber  disks,  where  depending  on  the  thickness  of  the 
disk,  they  either  observed  or  did  not  observe  cavitation,  and  two,  the  exper¬ 
iments  of  Ashby  et  al.(1989)  on  highly  constrained  lead  wires.  The  predictions 
of  the  model  have  been  consistent  with  the  observations.  In  Gent  and  Lindley's 
experiments  (see  description  previously)  they  studied  cavitation  in  rubber  disks 
of  various  thicknesses.  A  finite  element  stress  analysis  of  different  speci¬ 
mens,  as  the  load  is  increased,  yields  trajectories  in  stress  space  correspond¬ 
ing  to  the  evolution  of  the  stress  at  the  center  of  the  specimen,  as  shown  in 
Figure  13.  Cavitation  occurs  when  the  stress  trajectory  intersects  the  cavita¬ 
tion  curve.  The  figure  shows  that  the  two  thickest  specimens  do  not  intersect 
this  curve,  and  therefore  do  not  cavitate  at  the  center  of  the  specimen,  while 
the  others  do.  This  observation,  as  well  as  the  values  of  load  at  which  cavita¬ 
tion  does  occur  in  the  thinner  specimens,  agrees  with  Gent  and  Lindley’s  obser¬ 
vations  . 

Reference:  H-s.  Hou  and  R.  Abeyaratne,  Axi-symmetric  cavitation  in  solids,  manu¬ 
script  in  preparation. 
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Figure  1.  Variation  of  critical  stress  rcr  with  hardening  exponent  n  at 
fixed  Poisson’s  ratio  t/-0.3. 


Figure  2.  Variation  of  critical  stress  rcr  with  Poisson's  ratio  v  at 
fixed  hardening  exponent  n-2 . 
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Figure  6.  Pre-  and  post-bifurcation  sketches  of  a  material  element  in  a  cavita- 
ting  sphere. 
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Figure  9.  Cavitation  curve  0(r^,r 2)”0  for  neo-Hookean  material. 


Figure  10.  Cavitation  curve  ^(f^,»'2)'“0  for  elastic-plastic  materials  with 
different  hardening. 
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Figure  13.  The  Intersection  of  stress  trajectories  and  the  cavitation  curve  for 
rubber  disks  of  various  thicknesses. 


